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Overview: 

Research  investigations  supported  under  this  grant  addressed  problems  of 
solving  PDE’s  using  orthonormal  bases  with  controled  localization  in  the  time- 
frequency  domain  (e.g.  wavelets).  We  have  developed  a  number  of  fast 
algorithms  and  started  investigations  in  several  new  areas.  Among  the  results  is 
the  Unequally  Spaced  Fast  Fourier  Transform  (USFFT)  algorithm  and  a  new 
approach  to  numerical  homogenization. 

Research  Activities: 

With  James  Keiser  we  have  developed  an  adaptive  approach  for  solving  non¬ 
linear  evolution  equations  in  wavelet  bases.  The  analysis  provides  a  fundamental 
understanding  of  how  to  develop  algorithms  of  this  kind.  This  work  is  part  of 
Kaiser’s  thesis. 

Using  the  Battle-Lemarie  scaling  function  we  have  constructed  a  fast  algorithm  to 
evaluate  the  Fourier  transform  of  generalized  functions.  In  particular,  we  have 
developed  an  algorithm  for  the  Unequally  Spaced  Fast  Fourier  Transform  and 
tested  its  performance  in  dimensions  one  and  two.  Also  we  have  a  new  algorithm 
to  evaluate  the  Fourier  Transform  of  an  image  taking  into  account  the 
discontinuities  within  the  image.  Our  approach  also  has  implications  for 
numerical  solutions  of  PDE’s  which  we  plan  to  investigate  further. 

Jointly  with  Mary  E.  Brewster,  we  have  constructed  a  numerical  scheme  for 
homogenization  of  linear  differential  and  difference  equations  via  multiresolution 
analysis.  This  work  has  potential  for  significant  impact  by  permitting  us  to 
construct  the  equations  for  computing  the  projection  on  a  sparse  scale  of  the 
solution  of  equations  given  on  the  very  fine  scales. 

Jointly  with  John  Dunn  and  David  Gines  we  developed  a  practical  hybrid  method 
for  solving  static  and  quasi-static  problems  in  electromagnetics  using  wavelets. 


Jointly  with  Bruno  Torresani  we  developed  a  method  of  applying  operators  in 
wavelet  bases.  This  method  is  different  from  both  standard  and  non-standard 
forms  that  have  been  developed  previously  and  is  very  efficient  for  convolution 
operators,  e.g.  the  Hilbert  transform. 

Work  with  graduate  students: 

This  work  was  supported  via  an  AASERT  Grant.  I  am  working  with  two  students, 
James  Keiser  and  Robert  Cramer. 

The  work  with  James  Keiser  on  solving  partial  differential  equations  in  wavelet 
bases  is  almost  finished  and  James  will  be  graduating  soon.  The  results  will  be 
submitted  for  publication. 

The  work  with  Robert  Cramer  on  the  multidimensional  algorithms  for  the  fast 
evaluation  of  a  class  of  integral  operators  on  functions  is  continuing.  We  expect 
to  release  some  of  the  results  shortly. 
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